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Abstract — An overview on the present observational status and 
phenomenological understanding of cosmic rays above lO^'' eV is 
given. Above these energies the cosmic ray flux is expected to be 
gradually dominated by an extra-galactic component. In order 
to investigate the nature of this transition, current experimental 
activities focus on the measurement of the cosmic ray flux and 
composition at the 'ankle' or 'dip' feature at several EeV. At 
the ultra high energy end of the spectrum, the flux suppression 
above 50 EeV is now well established by the measurements of 
HiRes and the Pierre Auger Observatory and we may enter the 
era of charged particle astronomy. 

1. Introduction 

The all particle spectrum of cosmic rays is known to follow 
a power law, dN/dE cc E ^ over many orders of magnitude. 
However, at the highest energies, shown in Fig. [l] it exhibits 
three remarkable features. The knee [1], a steepening of the 
flux by A7 « 0.5, at a few PeV followed by a flattening called 
the ankle [2] at several EeV and a flux suppression at ultra high 
energies [3]. 

The first two features are suspected to be an indication of 
the end of the galactic cosmic ray spectrum and the transition 
to an extra-galactic component. 

At energies above several hundreds of TeV the particle 
fluxes are too low to allow for a direct measurement of the 
properties of cosmic rays. Instead, as will be explained in 
Sec. nil the analysis of air showers plays a crucial role to 
measure their flux and composition. These two observables are 
essential to study the transition from galactic to extra-galactic 
cosmic rays and to distinguish between the various models put 
forward to explain the ankle (see Sec Hill and lIVIl. 

At the same hour this talk was given, the first proton beams 
were injected to the Large Hadron Collider [4], that will 
eventually be able to accelerate protons up to 7 • lO^'^ eV. 
The ultra-high energy frontier of physics is however beyond 
10^" eV, where several cosmic rays have already been detected 
[5]-[7]. At these extreme energies, particles are expected to 
suffer significant energy losses during their propagation to 
earth. The corresponding flux suppression was predicted over 
forty years ago [3] and it is only now, that experiments 
gathered enough statistics to study it carefully (see Sec. IVli. 

The astrophysical sources that are able to accelerate par- 
ticles to such tremendous energies are still unknown. Their 
unambiguous identification requires to study the arrival direc- 
tions of cosmic rays, i.e. to do particle astronomy. 

Whereas the experimental knowledge on cosmic rays made 
a major leap forward in the current hybrid era [8], new 
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Fig. 1. All particle flux of cosmic rays ( [10] and references therein) 

projects aim to accumulate more data below 1 EeV and at ultra 
high energies. Moreover, new data on hadronic interactions 
at man-made particle accelerators are needed to facilitate the 
interpretation of air shower data [9]. These efforts will be 
described in Sec. |Vll 

II. Air showers 

Cosmic particles entering the earth's atmosphere sooner or 
later collide with the nuclei of the ak and initiate a particle 
cascade, the so-called air shower Since the thickness of the 
atmosphere is more than 20 radiation and interaction lengths at 
vertical incident, it constitutes a suitable calorimeter to study 
the properties of the primary cosmic ray particles. 
Air shower detectors either measure the lateral densities of 
particles at ground or the longitudinal development of the 
cascade in the atmosphere. 

The qualitative relation of these experimental observables 
to the energy Eq and mass A of the primary particle can be 
easily understood within the simple Heitler-model [11]-[14] 
of air showers. Here one assumes that after each hadronic 
interaction length. A, 7r-mesons are produced with an average 
multiplicity of (n). In each interaction an energy fraction of 
/ « 1/3 goes to neutral pions which decay immediately into 



two photons and thus feed the electromagnetic component 
of an air shower, that develops through pair-production and 
bremsstrahlung until the energy of electromagnetic secondaries 
falls below the critical energy eom (~ 81 MeV in air). In that 
way, the fraction of the primary energy in the electromagnetic 
component at the nth interaction at depth n ■ A increases to 
1 — (1 — /)" until the energy of the charged pions falls below 
the critical energy ech at which their decay length becomes 
smaller than the interaction length. If one furthermore assumes 
that a nucleus with mass A and energy is equivalent to A 
nucleons of energy Eq/A (the so-called superposition model), 
the following important relations can be deduced from this 
simplistic model: 

The average depth at which the electromagnetic cascade 
reaches its maximum, (X,„ax)i grows logarithmically with the 
energy per nucleon: 

(Xmax) =a + blg [{Eo/e,„,)/A] , (1) 

where the constants a and b depend on the properties of 
hadronic interactions, /, A and (n). Since most of the energy 
of the primary particle eventually ends up in the electromag- 
netic cascade, the integral electron number is a good estimator 
for Eq. Hence, the observation of the longitudinal development 
of an air shower with for instance fluorescence detectors [15], 
allows to measure simultaneously the primary energy and 
^max and can therefore be used to determine the absolute 
value of the average nuclear mass as a function of energy (cf. 
SecHYll. The elongation rate [16]-[18], d{X^^^) /A([gE) ot 
—d{\gA)/d{\gE), can be used to study the change of the 
primary cosmic ray composition with energy. 

The number of muons, iV^, from the decay of charged pi- 
mesons that can be detected by particle detectors on ground 
is given by 

N^^{E^le,^y' A^-^ (2) 

where again the properties of hadronic interactions are hidden 
in a single number, /3, that is proportional to the logarithm of 
the charged meson multiplicity. Given a similar relation for 
the number of electrons on ground, a detector that is capable 
to distinguish muons and electrons can therefore disentangle 
the energy and mass of cosmic primaries on a statistical basis. 

The simple Heitler approach is very useful to understand the 
principles of air shower physics, but of course in practice ex- 
periments employ full Monte Carlo simulations of air showers 
with for instance CORSIKA [19] to interpret their data. These 
simulations are however of limited predictive power, as they 
need to rely on models of hadronic interactions at energies 
beyond man-made accelerators. The related uncertainties are 
the source of considerable systematic uncertainties for the 
interpretation of air shower data (see [20], [21]). 

Since the energy estimated from the integral of the electro- 
magnetic longitudinal air shower development depends only 
very little on details of hadronic interactions [22], [23], modern 
air shower arrays like the Pierre Auger Observatory [24] or 
Telescope Array [25] use a hybrid approach to calibrate the 
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Fig. 2. Energy calibration of the surface detector of the Pierre Auger 
Observatory [26]: Shower size, 533, as a function of energy, -Epo, measured 
with the fluorescence detector. 

energy scale of their surface detector with the energy estimate 
of a fluorescence detector. An example is shown in Fig. |2l 
where the expected power-law dependence of the number of 
ground particles on the primary energy (Eq. (|2]i) can be seen. 

III. From Galactic to Extra-galactic Cosmic Rays 

The standard explanation for the knee-feature in the cosmic 
ray all-particle flux between 10^^ and 10^^ eV is that it marks 
the beginning of the end of the galactic cosmic ray spectrum 
due to the escape of the high energy charged particles from 
the magnetic confinement within the galaxy and/or the reach 
of the maximum energy of galactic accelerators (presumably 
supernova remnants). Although the current experimental data 
can not rule out alternative explanations for the knee (eg. [27]- 
[31]), extrapolations of the low energy cosmic ray data with 
a rigidity dependent cut off <x Z ■ Ec can describe the 
existing data very well [32], [33]. Moreover, the deconvoluted 
galactic mass spectra measured with KASCADE [34] show 
distinct knees for each elemental component, compatible with 
a rigidity dependent knee. 

Since Ec is of the order of PeV, it follows that above energies 
of several 10^* eV, the detected cosmic particles must be of 
extra-galactic origin. As a corollary, this assumption explains 
the lack of an observation of a strong anisotropy that would 
be expected for charged particles with a large gyro-radius at 
this energy. 

In the following we will describe three different models 
of the transition from galactic to extra-galactic cosmic rays 
(see [35]-[38] for recent reviews on this topic). Common to 
all of these models is that the properties of extra-galactic 



cosmic ray sources are described by just four parameters: the 
source emissivity (needed to adjust the overall normalization), 
the spectral index 70 of the energy spectrum at the source, 
the maximum energy Z ■ the source is able to accelerate 
particles to and the cosmological source evolution parameter 
m, that describes the source density n as a function of 
redshift z, dn/dz cc (1 + z)™. In order to simplify the 
discussion, we will restrict ourselves to the uniform source 
distribution model (m = 0) and assume that Ec is large 
enough to have no observational consequence within the 
statistical precision of current experiments. Furthermore we 
will assume the simple phenomenological rigidity-dependent 
parameterization from [32] to describe the 'standard' galactic 
cosmic ray component. Since there is an obvious disagreement 
of the measured ultra-high energy spectra [39]-[41] (cf. Fig[T]i, 
their energy scale will be 'adjusted' accordingly within the 
quoted systematic uncertainties. 

A. Dip Model 

In the so-called dip model [42], all extragalactic cosmic 
rays are assumed to be protons (at least after they escaped 
from their sources [43]) and the ankle feature is caused by 
energy losses suffered during the propagation to earth. 

First of all, the expansion of the universe causes adiabatic 
losses that are important for very distant sources (cf. Fig. [3^). 
Since it affects all energies equally, it does not change the 
spectral index observed at earth. 

A more important effect is the interaction of the cosmic ray 
protons with the photons of the cosmic microwave background 
radiation. At lower energies, the production of e+c^ -pairs 
through the Bethe-Heitler process is the dominant source of 
energy loss (cf. Fig. ^) and at energies above 10^^^ eV, 
the photon-proton center of mass energy is large enough for 
resonant photo-pion production, that gives rise to large energy 
losses [3] even for very close source (cf. Fig. [3};). 

As can be seen in Fig. |4^, the dip model can describe the 
data rather well. The transition between the galactic and extra- 
galactic cosmic ray component is at low energies just above 
0.1 EeV and produces the somewhat less prominent feature 
called the second knee [44]. 

Thus the dip model is a very economic approach in terms of 
assumptions, as it can explain all features in the cosmic ray 
energy spectrum in terms of the well understood interactions 
of protons with photons and predicts a transition energy that 
is low enough to be compatible with current estimates of the 
maximum energy of galactic accelerators [45]. 
However, the hard injection spectrum at the source is prob- 
lematic in terms of the overall energy luminosity of extra- 
galactic sources if extrapolated to low energies. Therefore an 
'artificial' softening of the spectrum below a certain energy is 
usually introduced [46] (not shown in Fig. |4^). 

B. Mixed Composition Model 

The dip-model works only for a pure proton beam since 
an admixture of heavier nuclei with a fraction of > 15% 
diminishes the agreement with the data considerably. This is. 
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(b) redshift and pair-production 
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Fig. 3. Fractional energy at earth, (i?carth/£^o), of protons with initial 
energy Ea as a function of the source distance/redshift (CRPROPA [47] 
calculation). 



because the threshold for e+e^ -production is proportional to 
the energy per nucleon and thus only relevant for protons at 
the energies of the ankle. Instead, cosmic ray nuclei loose 
their energy predominantly due to photo-disintegration at the 
giant dipole resonance [48]. The mean free path for photo- 
disintegration scales with the Lorentz-factor of the particle and 
drops rapidly above F > 10^. 

In mixed composition models [49], [50] the extra-galactic 
cosmic ray composition is assumed to be equal to the one 



measured at low energies in our galaxy. Due to the Lorentz- 
factor dependence of the energy loss, the individual spectra of 
nuclei with mass A are subsequently suppressed at energies 
above > A ■ lO^^ eV. 

As in case of the dip-model, this ansatz gives a good descrip- 
tion of the existing data (cf. Fig. I^'), but with a much softer 
extra-galactic source spectrum. The transition from galactic 
to extra-galactic cosmic rays is at about a factor 10 higher 
energies close to 1 EeV and correspondingly, this model needs 
galactic sources with a higher maximum acceleration energy 
than the dip-model. 

C. Ankle Model 

Finally, the traditional way to reproduce the ankle-feature 
is to model it as the intersection of a flat extra-galactic com- 
ponent with a steep galactic component (see for instance [52], 
[53]). In that case, as can be seen in Fig.|4j;), the galactic cos- 
mic ray spectrum extends to energies well above several EeV 
and thus would require a significant modification of the simple 
rigidity model of the knee. 

IV. Composition of UHECRs 

All of the transition models explained in the last section give 
a similar good description of the measured cosmic ray spectra 
under very different astrophysical assumptions. Since they 
differ substantially in the predicted cosmic ray composition as 
a function of energy, this observable is the key to distinguish 
between the models. 

The mass composition estimated from surface detector 
observables is compatible with large contributions from heavy 
elements up to the highest energies (see [54], [55]). These 
estimates rely to a large extent on an accurate prediction of 
the number of muons (cf. Eq. |2]l in air showers. However, 
modern hadronic interaction models differ by as much as 30% 
in the number of muons on ground [56], [57]. Moreover, the 
application of air shower universality to data from the Pierre 
Auger Observatory suggests, that current air shower simula- 
tions systematically underestimate the number of muons [58], 
[59]. 

The maximum of the longitudinal development of the elec- 
tromagnetic component of air showers (cf. Eq. [U provides a 
composition sensitivity that is somewhat less dependent on the 
details of hadronic interactions. As can be seen in Fig. |5^, all 
current air shower models give similar predictions of (X^ax) 
between 0.1 and 10 EeV. It is however worthwhile noting 
that this might be a mere coincidence, since also the predic- 
tions of the depth of the shower maximum can be changed 
significantly, if some more drastic (though experimentally not 
excluded) modifications of the hadronic interactions at high 
energies are assumed [60]-[63]. 

The shower maximum can be directly measured by fluores- 
cence detectors, that can infer the longitudinal shower devel- 
opment from the observation of fluorescence and Cherenkov 
light emitted by the shower as a function of height [75]. 
Observations of the lateral distribution of Cherenkov light at 
ground and its pulse shape are sensitive to X,„ax as well [76]. 
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Fig. 4. Models of the transition from galactic to extra-galactic cosmic rays 
vs. measurements of the all particle flux [34], [39]-[41], [51]. 




(a) Comparison to simulations [64] for proton and iron induced air showers (b) Comparison to transition models (calculations taken from [68], bands 
using different hadronic interaction models [57], [65]-[67] indicate the uncertainty from interaction models). 

Fig. 5. Measurements of (Xmax) from Cherenkov [69], [70] and fluorescence [71]-[74] detectors 



Measurements of the average of X„iax over almost three 
orders of magnitude in energy are shown in Fig. together 
with predictions from air shower simulations for proton and 
iron primaries. As can be seen, the data are indeed showing a 
trend from a heavy composition at low energies towards a light 
one at high energies, as would be expected from the transition 
models introduced in the last section. There are however 
systematic differences between the different experiments. The 
HiRes data, for instance, is compatible with a pure proton 
composition if compared to the QGSJET prediction, whereas 
the data from the Pierre Auger Observatory favors a mixed 
composition at all energies. 

A direct comparison of the data to the (X„iax) predicted 
by the dip-, ankle- and mixed-composition model is shown in 
Fig. ISj). Obviously, none of the three models gives a satisfac- 
tory description of the data, neither in shape nor the absolute 
(^max) value, but note that mixed-composition models have 
in principle enough parameters to be adjusted to the data. 

Until now we only discussed the average value of the 
shower maximum. The distribution of Xmax can potentially 
constrain the mass composition of cosmic rays even better. In 
the naive superposition model, one would expect that nuclei 
with mass A have smaller shower-to-shower fluctuations 
by a factor of 1/^/A. Correctly accounting for nuclear 
fragmentation leads to somewhat larger fluctuations of 
nucleus-induced showers [77]-[79], but still the width of the 
-''^max distribution of iron showers is about a factor three 
smaller than that for proton (about 20 and 60 g/cm^ at 1 EeV 
respectively). The analysis of the X,nax distribution requires 
however a good understanding of the detector resolution 
and corresponding composition estimates from the Xmax 
fluctuations are still contradictory (for instance pure proton 
in [80] and mixed in [81] above 1 EeV). 



V. The End of the Cosmic Ray Spectrum 

More than forty years after Greisen, Zatsepin and 
Kuzmin [3] (GZK) predicted a suppression of the cosmic 
ray flux due to interactions with the cosmic microwave back- 
ground (CMB) radiation and its existence has now finally been 
established with high significance by HiRes and the Pierre 
Auger Observatory [26], [40] (cf. Fig|6]l. Furthermore, Auger 
reported an anisotropy of the arrival direction of cosmic rays 
above 60 EeV [82] and set a limit of < 2% on the fraction 
of photons above 10 EeV. The latter excludes most of the 
top-down scenarios, i.e. cosmic ray production in decays of 
ultra-massive particles (see eg. [83]), that were motivated by 
the absence of a GZK feature in the AGASA spectrum [39]. 
The onset of the anisotropy at about the same energy as the 
GZK threshold suggests that the suppression is indeed due to 
propagation effects and not because the maximum energy of 
the sources is reached: If sources are isotropically distributed 
on large scales, local anisotropics can not be detected in a 
transparent universe, but only if propagation losses limit the 
distance from which cosmic rays can reach earth (the so-called 
GZK-horizon). 

It is a curiosity, that the thresholds for photo-pion production 
of protons with photons of the CMB is at a similar energy as 
the giant dipole resonance for iron nuclei. As can be seen in 
Fig. |6] the current statistical precision of the flux measure- 
ments at ultra-high energies is not sufficient to distinguish 
between the predictions for the spectral shape for cosmic 
rays with a pure proton [80] and iron [84] composition at 
the source. (It is worthwhile noting that the measured spectra 
are not corrected for the corresponding experimental energy 
resolution and if a deconvolution was applied 'true' shape 
of the flux suppression would get steeper). A possible way 
to resolve this degeneracy in the near future would be the 
detection of photons [85] or neutrinos [86] originating from the 
decay of pions produced during the proton propagation (nuclei 
are expected to produce much less neutrinos [87]-[89]). 
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Fig. 6. Measured flux at ultra-high energies compared to predictions for 
propagated proton and iron primaries (lines adapted from [80] and [84]). 




Fig. 7. Arrival directions of cosmic rays (E > 60 EeV) measured by the 
Pierre Auger Observatory (open circles) within its acceptance (shaded area) 
and location of active galactic nuclei (small dots) [82]. 



The anisotropy reported by the Pierre Auger Observatory 
was established by correlation of the arrival direction of cos- 
mic rays and the location of active galactic nuclei (AGNs) [90] 
(cf. Fig|7| within 3.1°. This angular scale is compatible with 
the deflections expected for protons in the galactic magnetic 
field (see for instance [91]), but from the current statistics 
it is not possible to distinguish if the AGNs from [90] are 
indeed the source of ultra-high energy cosmic rays or just a 
tracer of the true source distribution like the super-galactic 
plane [92], the large scale structure of nearby matter [93] or 
even a few sources producing nuclei of intermediate mass [94] 
that are spread out by magnetic fields. Note that a follow-up 
analysis of HiRes did not show a correlation [95], which may, 
however, be explained by the different energy scales of the 
two experiments. 

VI. Outlook 

The last years have brought a wealth of new precise data on 
cosmic rays above the knee, especially at ultra-high energies 
collected by the southern part of the Pierre Auger Observatory. 
Since its construction was just finished this year, one can soon 
expect updated results with increased statistics. Its northern 
part is planned to be built in Colorado, USA, and will increase 
the exposure of the observatory by a factor of seven and 
provide full sky coverage for particle astronomy [96]. The 
space-borne experiments TUS [97] and JEIVI-EUSO [98] will 
observe air showers from space and investigate the region 
above the GZK cutoff. 

At intermediate energies, the hybrid Telescope Array [25] 
started data taking and will study the region around the ankle 
with fluorescence detectors and a scintillator array. Both, the 



Pierre Auger Observatory and Telescope Array aim at cov- 
ering the transition region of galactic to extragalactic cosmic 
rays down to 10^' eV by using fluorescence telescopes with 
enlarged field of views and shielded particle counters [99]- 
[101]. The construction of the low energy enhancement of 
the southern Auger site is almost finished and 'first light' 
is expected in early 2009. Finally, the end of the galactic 
cosmic ray spectrum is currently observed by IceTop [102] 
and KASCADE-Grande [103] down to energies of 10^^ eV. 

This new cosmic ray detectors will thus cover more than 
four orders of magnitude from the knee up to beyond the 
GZK cutoff. A number of laboratory experiments will provide 
additional measurements to lower the systematic uncertainties 
of the cosmic ray measurements: Many of the modern 
detectors use fluorescence detectors to calibrate their energy 
scale and the absolute value of the fluorescence yield in air 
is one of the major contributions to the current energy scale 
uncertainties of HiRes and Auger. It is currently re-measured 
under various atmospheric conditions by several groups [104]. 
Furthermore, in order to diminish the uncertainties of the 
hadronic interaction models employed to interpret the cosmic 
ray data, more data from controlled interactions at accelerators 
are collected. The NA61 experiment [105] at the Super Proton 
Synchrotron at CERN will measure pion-carbon interactions 
above 300 GeV that are important for the last stages of the air 
shower development [106] and LHCf [107], TOTEIVI [108] 
and CASTOR [109] at the Large Hadron Collider will 
provide data on particle production in the forward region and 
the proton-proton cross section at center of mass energies 
corresponding to 10^^ eV in terms of primary cosmic ray 
energies. 
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